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The complete knowledge of Laplacian eigenvalues and eigenvectors of complex networks plays an 
outstanding role in understanding various dynamical processes running on them; however, deter- 
mining analytically Laplacian eigenvalues and eigenvectors is a theoretical challenge. In this paper, 
we study the Laplacian spectra and their corresponding eigenvectors of a class of deterministically 
growing treelike networks. The two interesting quantities are determined through the recurrence 
relations derived from the structure of the networks. Beginning from the rigorous relations one can 
' obtain the complete eigenvalues and eigenvectors for the networks of arbitrary size. The analyti- 

cs^ i cal method opens the way to analytically compute the eigenvalues and eigenvectors of some other 

deterministic networks, making it possible to accurately calculate their spectral characteristics. 
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03 . I. INTRODUCTION 
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■ As an interdisciplinary subject, complex networks have received tremendous recent interest from the scientific 
community [l], 0, 13, 0, HI, because of their flexibility and generality in the description of natural and manmade 
systems. A central issue in the study of complex networks is to understand how their dynamical behaviors are 
influenced by the underlying topological structure [H, 0, 13] ■ I n various dynamical processes, the effect of network 
structure is encoded in the eigenvalues (spectra) and their corresponding eigenvectors of its Laplacian matrix. For 
instance, the synchronizability of a network is determined by the ratio of the maximum eigenvalue to the smallest 
nonzero one of its Laplacian matrix ff| 0] • Again, for example, for continuous-time quantum walks Q in a network, 
the quantum transition probabilities [9j between two nodes are closely related to the eigenvalues and orthonormalized 
eigenvectors of its Laplacian matrix, which also determine the resistance between a pair of nodes and the average 
resistance of all couples of nodes in a resistance network [I3,[llJ. Thus, the complete (exact) knowledge of Laplacian 
QO ' spectra and eigenvectors is very important for understanding the network dynamics. 

Recently, a lot of activities have been devoted to the study of the spectra of complex networks [TH, El, EH, EH , 
providing useful insight into the topological properties of and dynamical processes on networks. However, most 
previous related studies have been confined to approximate or numerical methods, the latter of which is prohibitively 
£C) • difficult for large networks because of the limit of time and memory. Moreover, notwithstanding its significance, 
' relevant research on eigenvectors of Laplacian matrix of complex networks is much less. 

In the present paper, we investigate the Laplacian eigenvalues and eigenvectors of a class of deterministic treelike 
networks, which are constructed iteratively [la ] . By applying the technique of graph theory and an algebraic iterative 
procedure, we derive recursive relations for the Laplacian eigenvalues and eigenvectors of the networks. The obtained 
recurrence relations allow one to determine explicitly the full Laplacian eigenvalues and eigenvectors of the considered 
networks of arbitrary iterations from those of its initial structure. 



II. MODEL FOR THE GROWING TREES 



Here we introduce a model for a class of deterministically growing trees (networks) defined in an iterative way [16| , 
which has attracted an amount of attention 0, E3 • We investigate this model because of its intrinsic interest and 
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FIG. 1: (Color online) Illustration of a deterministic uniform recursive tree for the special case of m — 2, showing the first 
several steps of growth process. 



its deterministic construction, which allows one to study analytically its Laplacian spectra and their corresponding 
eigenvectors. 

The deterministically growing trees, denoted by Ut it > 0) after t iterations, are constructed as follows. For t = 0, 
Uo is an edge connecting two nodes. For t > 1, Ut is obtained from Ut—i by attaching m (to is a positive integer) 
new nodes to each node in \Jt-\- Figure Q] illustrates the construction process of a particular network for the case of 
to = 2 for the first four generations. 

According to the network construction, one can see that at each step U (t% > 1) the number of newly introduced 
nodes is L(U) = 2 to (to + 1)** . From this result, we can easily compute the network order (i.e., the total number 
of nodes) N t at step t, 

t 

N t = L{U) = 2 (to + 1)* . (1) 

ti=0 

The considered networks have a degree distribution of exponential form. Their cumulative degree distribution 
-Pcum(fc), defined to be the probability that the degree is greater than or equal to k, decays exponentially with k as 
Pcum(k) — (to + l) _ ~sr 16]. Their average path length (APL), defined as the mean of shortest distance between all 
pairs of nodes, increases logarithmically with network order [l7j . Thus, the networks exhibit small- world behavior 

Notice that the particular case of to = 1 is in fact a deterministic version of the uniform recursive tree (URT) [l9j , 
which is a principal model [2(| [2l[ of random graphs [22| . As one of the most widely studied models, the URT is 
constructed as follows [l9j |: start with a single node, at each time step, we attach a new node to an existing node 
selected at random. It has found many important applications in various areas. For example, it has been suggested 
as models for the spread of epidemics [23| . the family trees of preserved copies of ancient or medieval texts [24[, 
chain letter and pyramid schemes [25| , to name but a few. The to = 1 case of the networks studied here has similar 
structural properties as the URT, thus, we call the considered networks expanded deterministic uniform recursive 
trees (EDURTs), which could shed light in better understanding the nature of the URT. 

After introducing the EDURTs, in what follows we will study the eigenvalues and their corresponding eigenvectors 
of the Laplacian matrices of the EDURTs. 
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III. LAPLACIAN SPECTRA AND THEIR CORRESPONDING EIGENVECTORS 

Generally, for an arbitrary graph, it is difficult to determine all eigenvalues and eigenvectors of its Laplacian matrix, 
but below we will show that for U t one can settle this problem. 

A. Eigenvalues 



As known in Eq. |T]), there are 2(m + 1)* vertices in U t - We denote by V t the vertex set of U t , i.e., V t — 
\yx, V2, ■ ■ ■ , «2(to+i)*}- Let A t = [dij] be the adjacency matrix of network Ut, where = ctji — 1 if nodes i 
and j are connected, fly = aji — otherwise, then the degree of vertex Vi is defined as d Vi = ^2 



Let 



dv 2{m+1)t ) represent the diagonal degree matrix of Ut, then the Laplacian matrix of Ut is 



D t = &\&g(d Vll d V2 , 
defined by L< = D* A*. 

We first study the Laplacian spectra of U t , while we leave the eigenvectors to the next subsection. By construction 
it is easy to find that the adjacency matrix A t and diagonal degree matrix D t satisfy the following relations: 
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where each block is a 2(m +1)* 1 x 2(m + 1)* 1 matrix and I is the identity matrix. Thus, according to the above 
expressions (for A t and D t ) and the definition of Laplacian matrix, we have the following recursive relation between 
L t and L t _i: 



L t = Dt-At 

/ L t _ x +ml -I -I 
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Then, the characteristic polynomial of L t is 

Pt(x) 
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where the elementary operations of matrix have been used. According to the results in [26|, we have 



P t (x) = (det ((a: - 1)1))™ x det x — m — 



x - 1 



I-L 



t-i 



Thus, Pt(x) can be recast recursively as follows: 

P t (x) = {x - l) 2m{m+1)t 



x P t _i{y{x)), 



(5) 



(6) 



(7) 



where y(a;) = x — m— -^r- This recursive relation given by Eq. |(7j) is very important, from which we will determine the 
complete Laplacian eigenvalues of Ut and their corresponding eigenvectors. Notice that P t -\{x) is a monic polynomial 
of degree 2(m + l) t_1 , then the exponent of in P t _i(y;(x)) is 2(m + l)* -1 , and hence the exponent of factor x — 1 
in P*(x) is 



2m(m + l) t_1 - 2(m + 1)* _1 = 2(m - l)(m + 1) 



t— l 



(8) 



Consequently, 1 is an eigenvalue of L t , and its multiplicity is 2(m — l)(m + 

Note that E/j has 2(m + 1)* Laplacian eigenvalues. We represent these 2(m + 1)* Laplacian eigenvalues as 
A*, A|, . . . , ^2(m+i)" respectively. For convenience, we presume A* < A 2 < • • • < ^%{m+i) 



and denote by E t the set of 



these Laplacian eigenvalues, i.e. E t — {X\, A|, . . . , A|, m+1 N t }. All the Laplacian eigenvalues in set E t can be divided 
into two parts. According to the above analysis, A = 1 is a Laplacian eigenvalue with multiplicity 2(m — l)(rn + l)* -1 , 
which gives a part of the eigenvalues of L t . We denote by E t the set of Laplacian eigenvalues 1 of Ut, i.e., 



E t -{1,1,1,. ..,1,1} 

V v ' 

2(m-l)(m+l)'- 1 



(9) 



It should be noted that here we neglect the distinctness of elements in the set. The remaining 4(m + l)*^ 1 
Laplacian eigenvalues of Ut are determined by the equation Pt-i(ip(x)) = 0. Let the 4(m + l) t_1 eigenvalues be 
A^, A 2 , . . . , ^4(„ l+ i)t-i ) respectively. For convenience, we presume A* < A 2 < • • • < A^ m+1 - lt _ 1 , and denote by P t * the 
set of these eigenvalues, i.e., P t * = {A*, X 2 , . . . , ^4( m+ i)t-i }• Therefore, the set of all Laplacian eigenvalues for U t can 
be expressed as E t = E t U P 4 * . 



According to Eq. 0, for an arbitrary element in Et—\, say A* 1 G E t -\, both solutions of x — m — -^j = A' - 



are 



in E*. In fact, equation x 



m 



A- 1 is equivalent to 



x 2 -(X\- 1 +m + l)x + Xl- 1 = 0. 



(10) 
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We use notations A- and A* +2 , to represent the two solutions of Eq. (fTOj) . since they provide a natural increasing 

order of the eigenvalues of Ut, which can be seen from the argument below. Solving this quadratic equation, its roots 
are obtained to be Af = f*i(A* _1 ) and A* +2 (»n+i)*-i = ^(A-" 1 ), where the functions ri(A) and ra(A) satisfy 

n(A) = i (A + m+ 1 - 7(ATto+ l) 2 - 4A) , (11) 

r 2 (A) = ~ (A + m + 1 + V(A + to + l) 2 - 4A) . (12) 

Substituting each Laplacian eigenvalue of Ut-i into Eqs. (fTTj) and (I12p . we can obtain the subset _E t * of Laplacian 
eigenvalues of Ut- Since i?o = {0,2}, by recursively applying the functions provided by Eqs. (fTTj) and (fl2|) . the 
Laplacian spectra of Ut can be determined completely. 

It is obvious that both ri(A) and ^(A) are monotonously increasing functions, and that they lie in intervals 

[0, 1) and (1, +oo), respectively. On the other hand, since ri(A) — 1 = i (\ + m — 1 — ^/ (X + to — l) 2 + 4rnj < 0, 

we have ri(A) < 1. Similarly, we can show that ^(A) > 1. Thus for arbitrary fixed A', ri(A) < 1 < ^(A') 
holds for all A. Then we have the following conclusion: If the set of Laplacian eigenvalues for Ut-i is E t -\ = 
{A^" 1 , A 2 _ , . . . , A^m+i)'- 1 }> ^cn solving Eqs. (fTTj) and (fT2"j) one can obtain the subset E% of Laplacian eigenvalues 

for U t to be E* = {A*, A|, . . . , A* (m+1)t _i}, where X\ < X\ < ■■■ < ~X t 2{m+1)t - 1 < 1 < A* (ro+1)1 _ I+1 < A 2(m+1)t _ 1+2 < 

• • • < A^ m+1 ^ t _i . Recall that E t consists of 2(to — l)(m + 1)* _1 elements, all of which are 1, so we can easily get the 

set of eigenvalue spectra for U t to be E t = E? U E t . 

From above arguments, it is easy to see that for the special case of rn = 1, all the 2* +1 Laplacian eigenvalues of 
Ut are fundamentally distinct, which is an interesting property and has never (to the best of our knowledge) been 
previously reported in other network models thus may have some far-reaching consequences. For other cases m > 1, 
some eigenvalues (e.g., 1) are multiple, which is obviously different from that of to = 1 case. 

It has been established that Laplacian eigenvalues have connections with many contexts in the theory of networks. 
For example, they are closely related to the number of spanning trees on complex networks (27j . It has been shown 
that the number of spanning tress on a connected network G with order N, N s t(G), concerns with all its nonzero 
Laplacian eigenvalues Xi (assuming Ai = and Xi ^ for i = 2, • ■ • , N), obeying the following expression [2j| 

1 N 

N at (G) = -[]A,. (13) 

Since Ut are trees for all parameter to, according to Eq. (|13p . the product of all nonzero Laplacian eigenvalues for Ut, 
denoted by At, should be equal to Nt, which can be confirmed from the following argument. For t = 0, by construction 
it is obvious that Ao = iVo = 2; for t > 1, according to Eq. (fl0|) . we can easily obtain the following recursive relation 
A t = (to + 1) A t _i, which combining with the initial value Ao — 2 leads to A t = 2(to + 1)' = Nt- This proves that 
our computation on the Laplacian eigenvalues for Ut is right. 



B. Eigenvectors 



Similar to the eigenvalues, the eigenvectors of L t follow directly from those of L t _i. Assume that A is an arbitrary 
Laplacian eigenvalue of Ut, whose corresponding eigenvector is v G R 2(m+1)t , where R 2 ( m+1 )* represents the 2(m + l)*- 
dimensional vector space. Then we can solve equation (A I — L t )t> = to find the eigenvector v. We distinguish two 
cases: A G E\ and A G E t , which will be separately addressed in detail as follows. 

For the first case A G E^ , we can rewrite the equation (AI — L 4 )t> = as 



/ (A - to)I - L t -! I I 
I (A -1)1 

I (A -1)1 
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Vl 
V 2 
V 3 
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= 0, 



V I ••■ (A -1)1/ V v m+1 ) 

where vector Vi (1 < i < to + 1) are components of v. Equation ()14|1 results in the following equations: 

((A - m)I t -i - L t -i)v! + v 2 H h v m+1 = 0, 

ui + (A-l)^ = (2<i<m + l). 



(14) 



(15) 
(16) 
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Resolving Eq. (|T6|) , we find that 



v, = — 



A- 1 



Substituting Eq. (JT7J) into Eq. (JT3J) we have 

A — m 



A- 1 



vi (2 <i <m + l). 



vi = 0, 



(17) 



(18) 



which indicates that V\ is the solution of Eq. (fTS)) while (2 < i < m + 1) are uniquely decided by V\ via Eq. (|17p . 
In Eq. ([71), it is clear that if A is an eigenvalue of Laplacian matrix L t , then /(A) = A — m 



A-l 



must be one 



eigenvalue of L t _i. [Recall that if A = A- e E%, then y(A*) = A* 1 for i < 2(m +1)' 1 , or </?(Af) = Aj_ 2 (m+i)*- 1 for 
z > 2(m + l)* -1 -] Thus, Eq. (fT8|) together with Eq. ((Tj) shows that i>i is an eigenvector of matrix Lt_i corresponding 
to the eigenvalue A — m — j^j determined by A, while 



V 3 

V v m+ x j 



( 



Vl 



(19) 



is an eigenvector of L t corresponding to the eigenvalue A. 

Since for the initial graph Uq, its Laplacian matrix Lo has two eigenvalues and 2 with respective eigenvectors 
(1, 1) T and (1, — 1) T ; by recursively applying the above process, we can obtain all the eigenvectors corresponding to 
A e E*. 

For the second case of A S E t , where all A = 1, the equation (A I — ~L t )v — can be recast as 



/ (1 - m)I - L t _i I I I 
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I 



V 



Vl 

v?, 



= 0, 



(20) 



V v m+ x J 



I 0/ 

where vector Vi (1 < i < m + 1) are components of v. Equation (|2H)l leads to the following equations: 

vi = 0, 



V 2 + V 3 



i=0. 



(21) 
(22) 



In Eq. (|2ip . V\ is a zero vector, and we denote by Vij the jth component of the column vector Vi. On the other 
hand, Eq. (|22| gives us the following equations: 



V2,l 
V 2 ,2 



V3,l 
V3,2 



V 



2,2(m+l)* 



V 



3,2(m+l)* 



V 



v m +i,i = 

Vm+1,2 = 



m+l^m+l)'- 1 



The set of all solutions to any of the above equations consists of vectors that can be written as 



/ V2j \ 

\ V m +ld / 



f- 1 \ 
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V o J 
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(23) 
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where k% t j , fe.j , ■ ■ ■ , fem-lj are arbitrary real numbers. In Eq. (j23|) . the solutions for all the vectors Vi (2 < i < m+1) 
can be rewritten as 



\ 



V "In / 



(-1 -1 

1 
1 
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1 J 
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3.1 



fcl,2 
^2,2 
^3,2 



^l,2(m+l) t " 1 ^ 
fe2,2(m+l) t_1 
^3,2(m+l)* _1 



Wn-l^fm+l)*- 1 / 



(24) 



where fcjj (l < i < m — 1; 1 < j < 2(m + l)*^ 1 ) are arbitrary real numbers. According to Eq. (|24p . we can obtain 
the eigenvector u corresponding to the eigenvalue 1. Moreover, it is easy to see that the dimension of the eigenspace 
of matrix L t associated with eigenvalue 1 is 2(m — l)(m + 

In this way, all eigenvalues and their corresponding eigenvectors of Ut have been completely determined in a recursive 
way. 



IV. CONCLUSIONS 



In this paper, we have investigated the Laplacian eigenvalues and their corresponding eigenvectors of a family of 
deterministically growing treelike networks that exhibit small-world behavior. Making use of the methods of linear 
algebra and graph theory, we have fully characterized the Laplacian eigenvalues and eigenvectors of the networks, all 
of which are recursively determined from those for the initial network. Interestingly, we showed that for a particular 
case [m — 1) of the networks under consideration, all its Laplacian eigenvalues are disparate. We expect our results 
to be interesting in some fields of networks, such as random and quantum walks on networks, the computation of 
the resistance between two arbitrary nodes in a resistor network, the dynamics of coupled oscillators on networks, 
and so on. We also expect that the computing methods of eigenvalues and eigenvectors used here might be extended 
to other types deterministic networks, e.g., deterministic small-world networks [2^. l30l] and deterministic scale-free 
networks [3ll. l32l. l33l. I34L l35l I36T] . 
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